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By Molecular Dynamics simulation of a coarse-grained bead-spring type model for a 
cylindrical molecular brush with a backbone chain of A*";, effective monomers to which 
with grafting density a side chains with effective monomers are tethered, several 
characteristic length scales are studied for variable solvent quality. Side chain lengths 
are in the range 5 < < 40, backbone chain lengths are in the range 50 < A^^f, < 200, 
and we perform a comparison to results for the bond fluctuation model on the simple 
cubic lattice (for which much longer chains are accessible, A^ < 1027, and which 
corresponds to an athermal, very good, solvent). We obtain linear dimensions of side 
chains and the backbone chain and discuss their A^-dependence in terms of power 
laws and the associated effective exponents. We show that even at the Theta point 
the side chains are considerably stretched, their linear dimension depending on the 
solvent quality only weakly. Effective persistence lengths are extracted both from 
the orientational correlations and from the backbone end-to-end distance; it is shown 
that different measures of the persistence length (which would all agree for Gaussian 
chains) are not mutually consistent with each other, and depend distinctly both on 
Nb and the solvent quality. A brief discussion of pertinent experiments is given. 
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I. INTRODUCTION 



Macromolecules which consist of a "backbone" polymer, to which flexible or stiff side 
chains are grafted, the so-called "bottle brush polymers" , find very great interest recently (see 
Refsi"- for reviews). Varying the chemical nature of both backbone chain and side chains, 
their chain lengths {Nb, N) and the grafting density cr, the structure of these cylindrical 
molecular brushes can be widely varied. I.e., their local "thickness" as measured by the 
cross-sectional radius Rc or linear dimensions of individual side chains can be varied as well 
as their local "stiffness" , traditionally measured by "the" persistence length Ijt—, and their 
effective contour length Lc- We here use quotation marks with respect to "the" persistence 
length, because there is evidence, at least for the case of very good solvent conditions, 
that a unique persistence length measuring the "intrinsic" stiffness of a polymer cannot 
be defined in the standard fashion^^iii. Now an intriguing observation^^ is the sensitivity 
of the large-scale structure of these bottle brush polymers to solvent quality: one finds a 
thermally induced collapse of single macromolecules from cylindrical brushes to spheres, in a 
very small temperature range, and it is speculated that these bottle brushes could be useful 
as building blocks of "soft nanomachines"— . We also note that biopolymers with comb-like 
architecture are ubiquitous in nature (such as proteoglycans^^, or the aggrecane molecules 
that play a role in the soft lubricating layers in human joints^, etc.), and probably in this 
context temperature and/or solvent quality (or pH value) are relevant parameters as well. 

In view of these facts, a comprehensive clarification of how the properties of bottle brush 
polymers depend on solvent quality clearly would be interesting. Although there are oc- 
casional experimental reports, how particular linear dimensions of these polymers scale in 
various solvents (e.g.—"—), we are not aware of a systematic study of this problem. While 
work based on self-consistent field theory (SCFT) predicted already very early on^^ that the 
side chain gyration radius Rgs scales as Rgs oc N'^^'^ for good solvents and Rgg oc N'^/^ for 
Theta solvents, there is now evidence from experiment, simulation and theory that these 
power laws apply if at all only for side chain length of the order of 10^, which are of 
no practical relevance: Experiments have only studied the range iV < 10^, and the range 
accessible in simulations^'^>^>22i-^ is similarly restricted. Numerical modeling applying the 
Scheutj ens-Fleer version of SCFT has given clear evidence^l that even within this mean- 
field approach ^ 10^ is needed to reach the regime where the predicted power laws^^^^ 
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apply. While the investigation of the scaling (self-similar) properties of bottle brushes with 
extremely long backbone chains and very long side chains may be a challenging theoretical 
problem, it is of little relevance for the experimentally accessible systems, and clearly not 
in the focus of the present paper. Both SCFT theoriesi^i^ and scaling theories22 consider 
ideal chains aiming to describe the behavior at the Theta temperature, while previous sim- 
ulations have almost exclusively considered the good solvent case only. Notable exceptions 
are studies of globule formation of bottle brushes with very short side chains {N = 4 — 12) 
under poor solvent conditions^ and the study of microphase separation of bottle brushes 
with straight backbones in poor solvent^i^^. 

The present work intends to make a contribution to close this gap, by presenting a 
simulation study of a coarse-grained bead-spring type model of bottle brush polymers where 
the solvent quality is varied from very good solvent conditions to the Theta point regime. In 
the next section, we shall describe the model and simulation technique, while in section III 
we present our numerical results. In Section IV, a summary is given, as well as an outlook 
on pertinent experiments. 

II. MODEL AND SIMULATION METHOD 

Extending our previous work on the simulation of bottle brush polymers with rigid 
backbones^i^, we describe both the backbone chain and the side chains by a bead-spring 
model^"— , where all beads interact with a truncated and shifted Lennard- Jones potential 
Uh:]{r) and nearest neighbors bonded together along a chain also experience the finitely 
extensible nonlinear elastic potential ?7fene(^), ^ being the distance between the beads. 



while Ui^j{r > Vc) = 0, and where Vc = 2.5(Tlj. The constant C is defined such that 
Ui^j{r = Tf.) is continuous at the cutoff. Henceforth units are chosen such that ^lj = 1, 
(Tlj = 1, the Boltzmann constant /c^ = 1, and in addition also the mass ttilj of beads is 
chosen to be unity. The potential Eq. ([T]) acts between any pair of beads, irrespective of 
whether they are bonded or not. For bonded beads additionally the potential ?7fene('") acts. 



Thus 




(1) 




< r < ro , 



(2) 
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while UpENEif > f'o) = oo, and hence tq is the maximal distance that bonded beads can 
take. We use the standard choice^^ ro = 1.5 and k = 30. Related models have been used 
with great success to study the glassification of polymer melts formed from short chains^ 
and to study the effects of solvent quality of polymer brushes on fiat planar substrates^. For 
such models of brushes on a planar substrate, the implicit solvent model (Eq. (1)) has been 
compared with models using explicit solvent molecules, and it was found that the results 
are very similar. 

Note that in our model there is no difference in interactions, irrespective of whether the 
considered beads are effective monomers of the backbone or of the side chains, implying 
that the polymer forming the backbone is either chemically identical to the polymers that 
are tethered as side chains to the backbone, or at least on coarse-grained length scales as 
considered here the backbone and side chain polymers are no longer distinct. There is also 
no difference between the bond linking the first monomer of a side chain to a monomer 
of the backbone and bonds between any other pairs of bonded monomers. Of course, our 
study does not address any effects due to a particular chemistry relating to the synthesis of 
these bottle brush polymers, but, as usually done^i^i^ we address universal features of the 
conformational properties of these macromolecules. 

There is one important distinction relating to our previous work— on bottle brush 
polymers with rigid backbones: following Grest and Murat^, there the backbone was taken 
as an infinitely thin straight line in continuous space, thus allowing arbitrary values of 
the distances between neighboring grafting sites, and hence the grafting density a could 
be continuously varied. For the present model, where we disregard any possible quenched 
disorder resulting from the grafting process, of course, the grafting density a is quantized: 
we denote here by a = 1 the case that every backbone monomer carries a side chain, a = 0.5 
means that every second backbone monomer has a side chain, etc. Chain lengths of side 
chains were chosen as = 5, 10, 20, and 40, while backbone chain lengths were chosen as 
Nb = 50, 100, and 200, respectively. 

It is obvious, of course, that for such short side chain lengths any interpretation of char- 
acteristic lengths in terms of power laws, such as Rc oc N"''^ , is a delicate matter, u^s being 
an "effective exponent" and characterizes only the specified range of rather small values of 
A^, and not the limit N ^ oo considered by most theories^'^^'^"— . Thus, the actual value 
of t'eff is of limited interest, it only gives an indication to which part of an extended crossover 
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region the data belong. However, we emphasize that (i) our range of nicely corresponds 
to the range available in experimentsi^'^^^—^'^'—. (ii) The analysis in terms of power laws 
with effective exponents is a standard practice of experimentalists in this context (e-g.-^i^). 

We recall that for linear chains the Theta temperature for the present (implicit solvent) 
model has been roughly estimated^ as Tg ^ 3.0 (note, however, that there is still some 
uncertainty about the precise value of Tg: for a similar model^ the correct value of Tg, 
Tg ^ 3.18 in this case, could only be established for chain lengths exceeding N = 200). 
Thus, in the present work we have thoroughly studied the temperature range 3.0 <T< 4.0. 
From previous work— on rather long chains in polymer brushes on fiat surfaces, using the 
same model ( Eqs. ([I]), ([2])) to describe the interactions, it is known that for T = 4.0 one finds 
a behavior characteristic for (moderately) good solvents. Very good solvent conditions could 
be obtained from a slightly different model that has extensively been studied for standard 
polymer brushes^^"^, where in Eq. ([1]) the cutoff is chosen to coincide with the minimum of 
the potential, Vc = (and then also T = 1 can be chosen for this essentially athermal 

model). Rather than carrying out simulations for bottle brushes using this model, we found 
it more appropriate to compare to the athermal version of the bond fluctuation model on 
the simple cubic lattice, which describes very good solvent conditions'"^ and has been used 
in our earlier work-*^"^. The use of this model has several advantages: (i) due to the 
fact that excluded volume constraints can be monitored via the occupancy of lattice sites, 
a very efficient implementation of the pivot algorithm in this Monte Carlo approach has 
become possible^. Therefore, very large bottle brush polymers can be equilibrated, up to 
= 1027, a task that is very difficult to achieve by Molecular Dynamics (MD) methods 
(ii) The extent to which the bond fluctuation model and MD results agree (for comparable 
choices of parameters) yields some insight to what extent the preasymptotic regime that 
we study is model dependent. Of course, a truly universal behavior (apart from amplitude 
prefactors) can only be expected for the asymptotic regime where the side chain length 

— )■ oo, that is not accessible in simulations or experiments. One could ask why we do not 
use versions of the bond fluctuation model where an effective attraction between monomers 
is included^ to study the effect of variable solvent quality in the framework of this model. 
The reason is that for temperatures slightly below the Theta temperature already practically 
frozen configurations of monomers occur, with several monomers next to each other blocking 
any possibility to move. Thus the convergence towards equilibrium then is extremely slow. 
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However, since the application of the athermal version of the bond fluctuation model 
(BFM) to the simulation of bottle brush polymers is well documented (Hsu and Paul^ have 
given a careful discussion of the effort needed for the BFM to sample equilibrium properties.) 
in the recent literature^i^i^, we do not give any details here. 

In the MD simulation, the positions rj(t) of the effective monomers with label i evolve 
in time t according to Newton's equation of motion, amended by the Langevin thermo- 
stai^--2^^^ 

r«Lj^ = -VU, mt)}) - 7§ + m , (3) 

where Ui {{rj{t)}) is the total potential acting on the z'th bead due to its interactions with 
the other beads at sites {rj{t)}, 7 is the friction coefficient, and Ti{t) the associated random 
force. The latter is related to 7 by the fluctuation-dissipation relation 

(f,(t) ■ f,(t)) = QkBT^6,,6{t - t') . (4) 

Following previous work—"— we choose 7 = 0.5, the MD time unit 

/ 2 \ 1/2 

being also unity, for our choice of units. Eq. ([3]) was integrated using the leap frog al- 
gorithnt^, with a time step 5t = O.OOGtlj, and utilizing the GROMACS package^. For 
the calculation of properties of the bottle brushes typically 500 statistically independent 
configurations are averaged over. 

Of course, for bottle brushes with large A*";, equilibration of the polymer conformations is 
a difficult problem. Since we expect that end-to-end distance and gyration radius Rg of 
the whole molecule belong to the slowest relaxing quantities, we studied the autocorrelation 
function of Rg to test for equilibration, 

{(Rl{t' + t)~Rl) (Rl{t')-Rl)) 

m = — z — ^^^^2 — ^ (6) 

{{Rlit')-Rl) ) 

Note that R^ means an average of R^{t') over the time t' in the MD trajectory (and the 
average (■ ■ ■) is computed by averaging over 500 statistically independent runs). Despite a 
substantial investment of computer time, 0(t) still exhibits significant fluctuations (remem- 
ber that quantities such as Re and Rg are known to exhibit a "lack of self- aver aging^i^" , 



irrespective of how large A^;, is). Fig. [Ugives some examples for (f){t). If the number of samples 
would be infinite and 6t — 0, we should expect a monotonous decay of (f){t) towards zero as t 
becomes large. Due to the fact that the number of samples is not extremely large, and 6t not 
extremely small, Fig. [T] gives clear evidence for noise that is still correlated. We see that for 
small the noise amplitude starts out at ±0.1, and the time scale on which the fluctuations 
of change sign is at t ISOOtlj in (a), t ^ IOOOOtlj in (b), and t ^ 20000rLj in (c). 
While in cases (a)(b) the (statistically meaningful) initial decay of (f){t) occurs so fast that it 
can only be seen on a magnified abscissa scale (inserts), we see that in (c) the initial decay is 
also much slower, and the associated time scale is of the same order as the time over which 
fluctuations are correlated. We have carefully considered 0(t) for all cases studied, and we 
have concluded that for our largest system studied (shown in Fig. [U^c), with a total number 
of 8 000 effective monomers) the statistical effort is not yet sufficient to allow meaningful 
conclusions on the overall linear dimensions of the bottle brush, while in all other cases the 
effort was judged to be sufficient. The damped oscillatory character of the relaxation seen 
particularly in Fig. [H^c) could be a matter of concern; we attribute this relaxation behavior 
to particular slow collective motions (breathing type modes) of the chain. 

Fig. [2] shows a small selection of snapshot pictures of equilibrated bottle brush polymers. 
From these snapshot pictures it is already clear that the side chains cause a significant stiff- 
ening of the backbone, at least on a coarse-grained scale, and that bottle brushes where 
Nfy is not very much larger than look like wormlike chains. This conclusion corroborates 
pictures generated experimentally (by atomic force microscopy or electron microscopy tech- 
niques, e.g.-"^*^), but this observation should not mislead one to claim that the Kratky-Porod 
wormlike chain model^ often employed to analyze such micrographs provides an accurate 
description of bottle brush polymers, as we shall see. 

Of course, both the bead-spring model and the BFM are idealizations of realistic comb- 
branched polymers, where also bond-angle potentials and torsional potentials are present and 
contribute to the local chain stiffness. Thus it is gratifying to note that nevertheless measured 
structure factors of real bottle brush polymers^ can be mapped almost quantitatively to 
their simulated BFM counterparts^, if the lattice spacing is fixed at a length of a few 
Angstroms, and about 3 chemical monomers are mapped onto two effective monomers of the 
BFM. Residual minor discrepancies may to some extent be due to solvent quality effects^^. 
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III. SIDE-CHAIN AND BACKBONE LINEAR DIMENSIONS AND 
ATTEMPTS TO EXTRACT "THE" PERSISTENCE LENGTH OF BOTTLE 
BRUSH POLYMERS 

Fig. 13] presents log-log plots of the normalized mean square gyration radius of the side 
chains (Rg^) /{llN) as a function of side chain length (/;, is the bond length between 
successive monomers), comparing data for grafting densities a = 0.5 and a = 1.0, and 
two temperatures (T = 3.0 and T = 4.0, respectively). Three different backbone lengths 
{Nfy = 50, 100, and 200) are included, but the dependence of the data on backbone length 
is not visible on the scale of the graph. Such a dependence on A^";, would be expected due 
to effects near the chain ends of the backbone: in Ref.— it was shown that there is less 
stretching of the side chains in the direction perpendicular to the backbone near the chain 
ends of the backbone than in the central part of the backbone. However, this effect is almost 
completely compensated by the fact that side chains near the backbone chain ends are more 
strongly stretched in the direction of the backbone (side chains then are oriented away from 
the backbone chain ends). This fact has already been seen for bottle brushes with rigid 
backbone but free ends^. 

We have also checked that for the bead spring model average side chain linear dimensions 
for bottle brushes with rigid and flexible backbones are identical for the bond fluctuation 
model, at least in the accessible parameter regime. For the bond fluctuation model a similar 
equivalence has been found also for the radial density profile of the monomers^. Fig. H] shows 
then the temperature dependence of (-R^^) for the different choices of A^, demonstrating that 
the side chain extension is independent of A*";, and agrees well between rigid and flexible 
backbones. The small deviations observable between the rigid and the flexible backbone 
cases for large A^ are due to residual systematic errors of the simulations for the flexible 
backbone. The present work therefore suggests that the equivalence between bottle brushes 
with rigid and flexible backbones carries over to chains in variable solvent quality, down 
to the Theta point (but we caution the reader that this equivalence will break down for 
poor solvents, in the regime of intermediate grafting densities where for rigid backbones 
pearl-necklace structures occur— ^2^). 

The straight lines on the log-log plots in Fig. [3] illustrate the empirical power law 

{R%) (X AT^^- (7) 
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consistent with corresponding experiments (seei"^ and references quoted therein) and we 
find that u^s decreases with decreasing solvent quahty and with decreasing grafting density. 
In the good solvent regime, u^s > v = 0.588^, the established value of the exponent 
describing the swelling of linear polymers under good solvent conditions, Rg{N) oc N" for 

— oo^>^. As expected, u^s is still distinctly smaller than the value predicted from scaling^, 
Ucs = 2z//(l + u) ^ 0.74 or SCFT— , 3/4. But the slight enhancement of Ucs with respect to 
its value for free chains [i^ = 1/2 at the Theta point, T = 3.0; u = 0.588 for good solvents) 
is evidence that the side chains interact with each other, which is a prerequisite for the 
expected induced stiffening of the backbone chain. 

We now turn to the backbone linear dimensions. Here we first caution the reader that in 
good solvent conditions we clearly also expect that the mean square end-to-end distance of 
a bottle-brush polymer satisfies the standard scaling relation 

{Rlf^) oc N^"" , ly = 0.588 , iVf, ^ oo , (8) 

but the longer the side chain length the larger Nb must be chosen such that Eq. ([8]) can 
be verified. We start the discussion with data from the bond fiuctuation model, where 
(for < 24) data up to A'f, = 1027 are available^^. Fig. \5\ hence shows a log-log plot of 
(R'^b) / i'^^b-^b'^) versus A^. The factor 2 in the denominator is arbitrarily chosen, since for 
Gaussian chains {u = 1/2) the result would simply be the persistence length Ip in units of 
the bond length^. Thus, it was suggested^^ that in the excluded volume case one could 
introduce an effective persistence length Ip^n^N) via the definition 

lp,R{N) = {Rl)/{2hN^n, m^oo. (9) 

Then Fig. [5] can be interpreted as a plot of lpfi{N)/li, versus A^ in the region where it is 
basically independent of N^. The quantity Ip^^^N) defined in this way is compatible with an 
effective power law, 

lp,R{N) (X A^^ , (10) 

but the effective exponent clearly depends on N^., if A^^ is not chosen large enough (Fig. ED . 
Thus, Fig. [5] clearly shows that lp.B.{N) is not a quantity characterizing the intrinsic stiffness 
of a bottle-brush polymer. For large A";,, the exponent ( seems to saturate at a value close 
to 0.65 and 0.74 for a = 0.5 and a = 1.0, respectively, but for chain lengths A";, in the range 
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from Nb = 67 to A^^ = 131 it is only in the range from ( = 0.41 to ( = 0.55 and from 
C = 0.48 to ( = 0.67 for a = 0.5 and cr = 1.0, respectively. 

Thus, it is no surprise that for the bead-spring model choosing Nfy = 50, 100, or 200 
values of this effective exponent ( in a similar range are found (Fig. [7]). Of course, one must 
be aware that there is no strict one-to-one correspondence^ between the meaning of chain 
lengths Nf, and N in different models: actually there may be the need for conversion factors 
^^(BFM)^^^(BS) jy(BFM) /jy(BS) between the bead spring (BS) model and the bond fluc- 
tuation model (BFM). The same fact is true when we compare simulations to experiments; 
e.g., the data of Rathgeber et al.— could be mapped to the BFM2^ implying an equivalence 
between N^""'^'^ = 400 and N^™'^ = 259 and between A^('=^p) = 62 and N^™^ = 48, for 
instance. We also note that with decreasing temperature (R'^f,) decreases, irrespective of N 
(and also ( decreases). Again we remind the reader that this parameterization of our data 
in terms of effective exponents (also used in related experimental work— only can serve 
to indicate the place in an extended crossover region to which the data belong. 

When we analyze the gyration radius of the backbone as a function of backbone chain 
length for the BS model (Figs. [H]) we find that over a restricted range of Nb one can fit 
the data by effective exponents again, {R'^bV^^ oc A''^"*', with 0.55 < Vcs < 0.95, and it 
is seen that Ucs increases systematically with side chain length A^, and Ucs for T = 4.0 is 
larger than for T = 3.0, for the same choice of A^. Of course, this variation of the effective 
exponent is just a reflection of a gradual crossover from the rod-like regime ((i?^^/^ oc A^;,) 
to the self-avoiding walk-like behavior of swollen coils {{R^)^^"^ oc A'^^, with v ^ 0.588) in the 
good solvent regime (T = 4.0) or Gaussian-like coils ((-Rp^^^ oc A"^^^^) for the Theta point 
(T = 3.0), respectively. Similar plots as in Fig. [8]but for the BFM are shown in Fig. [9](a)(b). 
Due to the availability of equilibrated data for much longer backbone chain lengths Nh for 
the BFM, the plot of {Rl^)/{2llN^'') versus A^";, (Fig.[9](c)) indeed shows that the data settle 
down to a horizontal plateau implying that the asymptotic region where Eq. (Q can be 
applied indeed is reached. Instead of R^f^ in Eq. ([9]), we use i?^^ such that the estimate of 
the A^-dependent effective persistence length Ipj^g is given by the horizontal line in Figj^c). 
Actually, for fixed side chain length A^ the different choices of Nh can be collapsed on a 
master curve, when A";, is rescaled by Sbiob-, which can be interpreted as an estimate for the 
number of segments per persistence length. In our previous work-iii, we have shown that 
bottle brushes under very good solvent conditions can be described as self-avoiding walks 
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of effective "blobs", having a diameter that is just twice the cross- sectional radius of the 
bottle-brush, and containing a number Sbiob of backbone monomers. Note that hence Sbiob is 
not a fit parameter, but has been determined independently. This data collapse on a master 
curve is demonstrated for the BFM for the case a = 1, = 6, 12, 18, and 24 (Fig, IHI^d)). 

Since "the" persistence length of polymers in the standard textbooks^"- is traditionally 
defined from the decay of bond orientational correlations along the chains, we next focus 
on this quantity. Defining the bond vectors in terms of the monomer positions as 

= rj-|_i — fi, i = 1, . . . , Nh — 1, this bond orientational correlation is defined as 

{cos9{s)) = 1^'^———— {ai-tti+s) (11) 

Note that (c^) = If and hence (cos6'(0)) = 1, of course. Considering the limit A'';, — oo, and 
assuming Gaussian chain statistics, one obtains an exponential decay, since then (cos6'(s)) = 
(cos6'(l))'* = exp [s ln(cos6'(l))], and thus 

(cos^(s)) = exp [sh/lp] , = -ln(cos^(l))/4. (12) 

However, it has been shown by scaling and renormalization group arguments^ and verified 
by simulationsi^iii that in the good solvent case there actually occurs a power law decay 

(cos^(s)) oc s"'^, /3 = 2(1 - z/) ^ 0.824, s ^ oo (13) 

while for chains at the Theta point^ or in melts^ one has 

{cose{s)) (X s-^/\ s^oo. (14) 

As far as bond orientational correlations are concerned, Gaussian chain statistics hence is 
misleading for polymers, under all circumstances. However, Shirvanyants et al.— suggested 
that for semifiexible polymers (at the Theta point) one can use still Eq. f|T2|) but only for 
1 < s < s*, where s* oc Ip/h controls the crossover from the simple exponential decay, 
Eq. f|T2|) . to the power law, Eq. f|T^ . Indeed, for a simple self-avoiding walks (SAW) model 
on the simple cubic lattice, where chain stiffness was controlled by an energy cost Ef, when 
the chain makes a 90° kink on the lattice^, it was shown that this suggestion^ works 
qualitatively, also in the good solvent case, with Ip oc exp{eb/kBT) for Sh/kBT ^ 1, although 
the crossover between Eqs. ( IT2|) and ( !T3|) is not sharp but rather spread out over a decade 
in the variable s. 
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Motivated by this findingii, Figs. [101 [H] hence present a few examples where (cos6'(s)) 
is plotted vs. s on a semi-log plot, to test for a possible applicability of Eq. f lT2|) for not 
too large s. However, we find that in fact Eq. 0121) does NEVER hold for small s {s = 1, 
2, 3), unlike the semifiexible SAW model of Ref.— , rather there occurs a very fast decay 
of (cos6'(s)) following a strongly bent curve (only two successive values of s can always 
be fit to a straight line, of course, but there is never an extended regime where a straight 
line through the origin, (cos6'(s = 0)) = 1, would be compatible with the data). It is 
also remarkable that this initial behavior is almost independent of a and A^; an obvious 
interpretation is that the stiffening of the backbone caused by the presence of long side 
chains is effective only on mesoscopic length scales along the backbone, but not on the scale 
of a few subsequent backbone bonds, which maintain a high local flexibility. Only for s > 4 
the data are compatible with a relation 

(cos6'(s)) = aexp(— 6s) , 4 < s < Smax , (15) 

where a and b are phenomenological constants, and s^ax depends on both A^, cr, and T 
distinctly (but cannot be accurately obtained from our simulations, because for (cos^(s)) < 
0.03 the statistical accuracy of the data deteriorates.) Obviously, the relation (Eq. ( IT2l) ) 
= — ln(cos 6{l))/£i, fails, but it seems tempting to identify an effective persistence length 
Ip^ as Ip^/ib = b^^- However, when we would define the persistence length in this way, we 
obtain the result that 1^^ depends on both A^e, and on T, not only on the side chain length A^ 
and grafting density cr (Fig. [12] and Table [Tj). For the bead-spring model, one often estimates 
that the length unit (ctlj = 1) physically corresponds to about 0.5 nm. The data for i'^{T) 
at cr = 1 and good solvent conditions for A^f, = 200 then would span a range from about 5 
to about 40 nm, i.e. a similar range as proposed in recent experiments^^. But already our 
previous work on the athermal bond fluctuation mode^>^^^ has given some evidence, that 
defining a persistence length from a fit of the data for (cos 6'(s)) to Eq. f[T5|) is not suitable 
to obtain a measure of the local intrinsic stiffness, since Ip^ depends on Ni,. The present 
data show that Ip^ depends on T as well, and hence is not just controlled by the chemical 
architecture of the bottle brush (via the parameters, a, N, and A^f,), but depends on solvent 
conditions as well. 

Thus, we argue that the physical significance of a persistence length 1°^ extracted from 
bond orientational correlations in this way is very doubtful, even at the Theta point. We 
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also note that sometimes, due to curvature on the semi- log plot (e.g. Fig. [TTT d)) such fits 
are ill-defined. 

Finally, we consider the possible validity of the Kratky-Porod result for the end-to-end 
distance of the chains at the Theta temperature, where for ^ oo we have (apart from 
logarithmic corrections) (Rl) oc Nf,, i.e. for this property a formula analogous to Gaussian 
chains holds. The Kratky-Porod result describes the crossover from rods to Gaussian chains, 

= 1 - I [1 - exp(-L//,)] . (16) 

We can define an effective exponent UcslL/lp) in terms of the logarithmic derivative of this 
function. Fig. [13] (a), (b). 

Using the data for u^s at Nf, = 100 from Fig. [8](a)(c) (and similar data not shown for (-Rg^) 
) we hence obtain estimates for the ratios L/lp for the various choices of for both a = 0.5 
and a = 1. Note that the errors in our estimation of u^s translate into rather large errors in 
our estimates for L/lp (both these errors of our data are indicated in Fig.[T3](a)(b)). Since we 
expect that the actual variation of (Rli) with A'^^ exhibits slight curvature (although this is 
hardly detected in Fig. [Hta)(c)), we use only this intermediate value of Ni, = 100 to estimate 
the relation between and L/lp, and we use neither Nf, = 50 nor A'^^ = 200 for this purpose. 

If the contour length L of the Kratky-Porod model, Eq. f|T6l) . would simply be the "chem- 
ical" contour length L^h = N^lf,, where the bond length l^ connecting two neighboring 
monomers along the backbone is ~ 1, our results for L/lp would readily yield explicit 
results for Ip. However, using then Eq. ( fT6l) to compute (i?g^)/L^ yields an overestimation 
by about a factor of five for cr = 1. This discrepancy proves that L is significantly smaller 
than Lch- Since Eq. (17) implies that there is a unique correspondence between u^s and the 
ratio ip/L, the factor of five discrepancy for {R^) (which is proportional to the product of ip 
and L) means that both ip and L must be smaller by the same factor (at about v^)- Thus 
we define L = Nblf^ , where lf^{< 1) is interpreted as the average projection of a backbone 
bond on the direction of the coarse-grained contour of the wormlike chain. Thus we also 
can use the data shown in Fig. [T3r a)(b) to obtain explicit estimates for lp{N), taking into 
account L = Nfylf^ instead of L = Nf,lf,. Interestingly, the mechanisms leading to If^ < 1 
is also evident in the presence of the constant a < 1 in the fits of (cos6'(s)) = aexp(— 6s) 
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in Figs. Uni [H Using then the estimates if {a = 0.5) ^ 0.33 and if (cr = 1) ^ 0.45, the 
data for (i?g^)/(2/p(A^)L) at T = 3.0 are roughly compatible with Eq. f|T6|) . when we use 
L/lp{N) from Fig. [T^ a)(b) and take L = IfN^, as is shown in Fig. [T^ c)(d). In this way, 
we have defined a correction factor a,. = [if) ^ in Fig. [T31 which was assumed to depend 
neither on nor on Nf,. Note that while only Nf, = 100 was used in Fig. [T^ a)(b). data for 
Nb = 50, and 200 are included in Fig. [T3]^c)(d). The resulting values of lp{N) which hence by 
construction do not depend on Nf, (Table [T]) are considered to be the most reliable estimates 
for the considered model. However, we emphasize that Eq. ( IT6|) for bottle brushes is useful 
only if the solution is at Theta conditions, but not in the good solvent regime. Using the 
data for u^s at A^";, = 100 from Fig. |Hl^a)(c) (and similar data not shown for {Rl^} ) we hence 
obtain estimates for the ratios L/lp for the various choices of for both a = 0.5 and a = 1. 

IV. CONCLUSIONS 

In the present paper, a coarse-grained bead-spring model for bottle brush polymers was 
studied by Molecular Dynamics methods, varying both the chain length of the backbone 
and of the side chains (A^), for two values of the grafting density, under variable solvent 
conditions. The main emphasis of the present work was a study of the various characteristic 
lengths describing the conformation of the macromolecule, contrasting the behavior under 
Theta conditions with the behavior in the good solvent regime. Also a comparison with 
corresponding results for the bond fluctuation model has been performed; this athermal 
model represents very good solvent conditions, and while for corresponding parameters it 
yields results that are rather similar to the results for the bead spring model in the good 
solvent regime, it has the advantage that much larger values of A^;, can be studied. 

Among the quantities that have been studied are the end-to-end distances of side chains 
and of the backbone, as well as the respective gyration radii, but also mesoscopic lengths 
that are particularly popular when a description in terms of the Kratky-Porod wormlike 
chain model is attempted, namely "the" persistence length and "the" contour length of the 
wormlike chain. We have shown that for the range of rather short side chain lengths A^ that is 
only accessible, either in simulation or experiment, one often finds that a description in terms 
of power laws with effective exponents is adequate. Specifically, the side chain radii vary as 
^ jyv^ff^ ^ Q gQ ^j^g rjnj^g^^ tempcrature T = Tg) or u^s ^ 0.63 to 0.66 (in 
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the good solvent regime), see Fig. |3l These effective exponents (for the present models with 
flexible backbones) are the same as for otherwise equivalent models with stiff backbones, and 
also are practically independent of the backbone chain length. While the resulting effective 
exponents are systematically smaller than the asymptotic values predicted theoretically for 
the limit — oo, v ^ 0.74 (good solvents) or u = 2/3 (Theta solvents), respectively, they 
are comparable to the results of pertinent experiments: Zhang et al.->^ report Ucs ~ 0.67 for 
a bottle brush where the chemical nature of the main chain and the side chains is identical. 
For other systems exponents in the range 0.56 < z/cff < 0.67 are reportecU^, again values 
comparable to our findings. Some experimental studies also report that somewhat different 
results occur when different solvents are usedii^ii^ii^ii^, but a systematic variation of solvent 
quality to our knowledge has not yet been performed. Our results have shown that for 
T > Tg the side chain radii vary only very weakly with temperature (which in our model 
causes the change in solvent quality), see Fig. HI 

Also the backbone end-to-end distance shows a power law variation with the side chain 
length, (-Rgf,) oc A^'' (Fig. O - [7]), reflecting thus the systematic backbone stiffening with 
increasing side chain length. In this case, however, the effective exponent depends on the 
backbone chain length A^, for A^^ < 300, before a saturation at C ~ 0.65 {a = 0.5) or 0.75 
(a = 1) occurs, in the case of the bond fluctuation model. Even somewhat larger exponents 
are observed for the bead-spring model (up to ( ^ 0.88, Fig. [7j). 

A further effective exponent describes the variation of the backbone end-to-end distance 
and gyration radius with backbone chain length (Figs. [H [9]). We have interpreted this 
variation as a crossover from rod-like behavior at small Ni, to self-avoiding walk-like behavior 
(under good solvent conditions) or random walk-like behavior ( under Theta conditions) with 
increasing Nh. 

Finally, "the" persistence length has been estimated in various ways, and demonstrating 
that the different estimates are not mutually consistent with each other it was argued that 
a unique persistence length does not exist. While a regime of exponential decay of the bond 
autocorrelation function often can be observed over some intermediate range of contour 
distances (Fig. [TDl [TT]) . the resulting estimates for a persistence length do not only depend 
on side chain length A^ and grafting density cr, but also on backbone chain length Ni, and 
on temperature T (Fig. [T2|) . However, when one studies the variation of the end-to-end 
distance of the backbone for T = Tg an analysis in terms of the Kratky-Porod wormlike 
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chain model {Eq. (fT6|) } becomes feasible (Fig. [T3|) . But one must not identify the contour 
length L implied by this model with the "chemical" contour length Lch = iVf,/;,, where Ih is 
the actual bond length, but rather one has L = Nijlf^ with If^ distinctly smaller than 1^. 
This effect results from the flexibility of the backbone on small scales; only on the scale of 
several backbone bonds does the stiffening due to the mutual side chain repulsions come 
into play. Thus, at the Theta point both an effective contour length L and a persistence 
length ip{N) are well-defined quantities, in terms of a fit of the data to the Kratky-Porod 
model, while under good solvent conditions such an analysis is not appropriate. 

We hope that the present work stimulates some more systematic experimental work on 
these issues. Of course, for real bottle brushes where the chemical nature of the backbone 
often differs from the side chains, the interesting possibility arises that the backbone and 
the side chains have rather different Theta temperatures. Another subtle issue concerns 
the precise location of Theta conditions. While we expect that for the limit A^";, — oo at 
fixed ratio N/Ni, the Theta conditions coincide with the Theta conditions for simple linear 
polymers {Ni, — )■ oo, = 0), it is not clear whether the location of the Theta point for the 
limit A^b — oo at fixed (small) A^ depends significantly on A^. Thus, also more theoretical 
work clearly is needed. 
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TABLE I. Estimates of the effective persistence length 1^^ = b~^ih shown in Fig. [T2]are hsted for 
the BS model with the grafting densities a = 0.5 and 1.0 at temperatures T = 3.0 and 4.0. Various 
values of the backbone length Ni, and the side chain length are chosen here. All length quoted 
in this table are given in units of £b (with Ih ~ 0.97cjlj). Note that £^ depends not only on N and 
cr, but also on Ni, and T, and hence does not seem as a characteristic of intrinsic chain stiffness. 
Estimates of lp{N) shown in Fig. [13] are also listed for comparison. 
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(a) (b) 




FIG. 1. Plot of 4>{t) (see Eq. ([6])) versus the time t in units of tlj. Several cases are shown: 
a = 0.5, Nb = 100, iV = 10, T = 3.0 (a) a = 1.0, Nb = 100, N = 20, T = 3.0 (b); and a = 1.0, 
Nb = 200, N = iO, T = 3.0 (c). The fluctuations of (j){t) change sign initially at t ^ ISOOtlj (a), 
t S3 IOOOOtlj, and t « 20000tlj (c), as shown in the inserts. 
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FIG. 2. Selected snapshot pictures of equilibrated configurations of bottle brush polymers. Back- 
bone monomers (when visible) are displayed in magenta color, side chain monomers in green. 
Cases shown refer to cr = 0.5, Nb = 100, N = 10, T = 3.0 (a) and T = 4.0 (b), as weU as a = 1.0, 
Nb = 100, iV = 20, T = 3.0 (c) and T = 4.0 (d). 
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FIG. 3. Log-log plot of mean square radius of gyration of the side chains normalized by their 
chain length and by the square bond length between successive monomers, {Fl^g) / {if^N), versus 
side chain length N , for a = 0.5 (a) and a = 1.0 (b). Data for the bead spring model at two 
temperatures (T = 3.0 and T = 4.0), as indicated in the figure are included, as well as three 
backbone chain lengths [Ni, = 50, 100, and 200, respectively). For comparison, also data for the 
athermal bond fluctuation model for comparable backbone chain length are included. Straight 
lines indicate effective exponents fcff ^ 0.60 (T = 3.0) or u^q ~ 0.63 (T = 4.0) in case (a), and 
~ 0.60 (T = 3.0) or u^s ~ 0.64 (T = 4.0) in case (b). For the bond fluctuation model under 
very good solvent conditions the slightly larger effective exponent {ucs w 0.65 for a = 0.5, and 
feflF ~ 0.66 for fj = 1.0) than for the off-lattice model with T = 4.0 results. 
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FIG. 4. Temperature dependence of the normalized mean square radius of gyration {R'^^) /{I'^N) 
for four chain lengths (A^ = 5, 10, 20, and 40), and three backbone lengths (A'^f, = 50, 100, and 200), 
respectively. All data are for the bead spring model and a = 1. Data for bottle-brush polymers 
with flexible backbone are connected by curves to guide the reader's eyes. Data are only shown for 
T = 3.0 and 4.0 for the rigid backbone case. 
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FIG. 5. Log-log plot of [W-^^ I {21\NI") = ip^n/ib versus side chain length N, for the bond fluctua- 
tion model with a = 0.5 (a) and a = 1.0 (b) including different choices of A'^^, as indicated. 
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FIG. 6. Plot of the effective exponent extracted from the slope of the straight lines in Fig. [U 
versus the backbone chain length Ni, for a = 0.5 (a) and a = 1.0. All data are for the athermal 
bond fluctuation model. 
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FIG. 8. Log-log plot of {Rl^^l'^lh versus backbone chain length A'^f,, for a = 0.5 (a)(b) and a = 1.0 
(c)(d), including both T = 3.0 (a)(c) and T = 4.0 (b)(d), for the bead-spring model. Several side 
chain lengths N are included as indicated. Straight lines indicate a fit with effective exponents i^eff- 
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FIG. 9. Log- log plot of for the bond fluctuation model versus backbone chain length 

Nil, for < < 259, and several side chain lengths N , for a = 0.5 (a) and a = 1.0 (b). Part (c) 
shows the plot of rescaled radius of gyration (^^5) / {'^I'b^b'^) °^ bottle brush polymers versus 
Nh for a = 1.0. The persistence length Ip^Rg is determined by the values of plateau. Part (d) 
shows a crossover scaling plot, collapsing for = 6, 12, 18, and 24, a = 1.0, and all data for 
67 < A'f, < 1027 on a master curve, that describes the crossover from rods (^{R'^h)^^'^ oc A'^;,^ to 
swollen coils (^(-Rg^)^/^ (x , u ^ 0.588^ . For this purpose. A*";, is rescaled with the blob diameter 
■Sbiobi which has been determined to be Sbiob = 6, 10, 12, and 14 for = 6, 12, 18, and 24, 
respectively^iii. 
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FIG. 10. Semi-log plot of (cos0(s)) vs. s for the bead spring model with a = 0.5. Data are chosen 
for Nb = 50 (a)(b) and Nb = 200 (c)(d) at T = 3.0 (a)(c) and T = 4.0 (b)(d), respectively. Several 
choices of N are shown as indicated. Straight lines illustrate fits to aexp(— 6s), with constants a, 
b quoted in the figure. 



(a) (b) 



A 




FIG. 11. Semi-log plot of (cos^(s)) vs. s for the bead spring model with a = 1.0 and = 100 
at the temperatures T = 3.0 (a) and T = 4.0 (b), respectively. Several choices of N are shown as 
indicated. Straight lines illustrate fits to a exp(— 6,s), with constants a, b as quoted in the figure. 
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FIG. 12. Log- log plot of the effective persistence length /it = as a function of side chain 
length N for T = 3.0 (a)(c) and T = 4.0 (b)(d) for several choices of Nb. Results for a = 1.0 and 
(T = 0.5 are shown as indicated and the exponent estimates C extracted from the slope are listed. 
Values of Ip^ are also listed in Table HI 
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FIG. 13. Plot of the function i'cf[{L/lp) versus L/lp as predicted from the Kratky-Porod model, 
Eqs. (jl6p . (jl7p (full curves). The numbers for L/lp extracted for a = 0.5 (a) and a = 1.0 (b) are 
quoted in the figure for various at the abscissa. The average of the exponents fcff from Fig.[8|^a)(c) 
and from R^^ (not shown) are quoted on the ordinate. Log-log plot of ar{R^fj) / {2lp{N)L) vs. 
L/lp{N) for a = 0.5 (c) and a = 1.0 (d). Several choices of are shown as indicated. Data 
plotted by the same symbol correspond to N^, = 50, 100, and 200 from left to right. The prediction 
for the Kratky-Porod model {Eq. (jl6p } is also shown for comparison. Approximate data collapse is 
obtained by introducing a factor (a^ = 9.2 it 1.4. and 4.9 it 0.2 for a = 0.5 and 1.0, respectively), 
cf. text. Values of lp{N) are also listed in Table HI 
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